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Abstract Given a parametrised weight function u(x,n) such that the quo- 
tients of its consecutive moments are Mobius maps, it is possible to express 
the underlying biorthogonal polynomials in a closed form ||. In the present 
paper we address ourselves to two related issues. Firstly, we demonstrate 
that, subject to additional assumptions, every such to obeys (in x) a linear 
differential equation whose solution is a generalized hypergeometric function. 
Secondly, using a generalization of standard divided differences, we present a 
new explicit representation of the underlying orthogonal polynomials. 
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1 Introduction 



Let us consider a one-parametric family of weight functions lu(x, fi) > 0, x > 0, fix < ft < 
We define a biorthogonal polynomial p n E P n x (^_,/z + ) n , p n ^ 0, by means of the 
orthogonality relations 
-l 

p n (x;(j,i,[j l 2,...,[j, n )u(x,[xt)dx = Q, £ = 1,2, ...,n, (1.1) 



where /Ui, fi 2 , ■ ■ ■ , fin ar e distinct. (Here - and elsewhere in this paper - P n denotes the set 
of all nth degree polynomials.) 

The theory of general biorthogonal polynomials has been described in detail in [Q] 
and further developments have been reported in |5|, g]. It possesses several interesting 
applications, ranging from numerical methods for ordinary differential equations || |8| to 
approximation theory p| to location of zeros of polynomial transformations and orthog- 
onal expansions [||, ||, |1 lfl . An exposition of biorthogonal functions and their applications 
features in Claude Brezinski's monograph Q. Herewith we revisit these elements of the 
theory that are germane to the work of the present paper. 

Biorthogonal polynomials exist and are unique (up to a nonzero multiplicative con- 
stant) if and only if 



det 



for every distinct fi\ , fi2 , 



miGui) 
tni(// 2 ) 



TOm-l(Ml) 

m m _i(// 2 ) 



) mi(/^ 
. ,/% G (fi~,fi + ), where 



^0 



x e uj(x, fi) dx 



is the £th moment of to, £ € Z + Q. The zeros of p n need not, in general, reside in (0, oo) 
However, all the zeros are in the support if to is strictly sign consistent (SSC), i.e. if 



det 



u(xi,n\) u(x 2 ,fix) ■■■ u(x m ,fix) 
u{xi,H2) uj(x 2 ,fi 2 ) ■■■ uj(x m ,fi 2 ) 



7^0 



u{xx,u. m ) u(x 2 ,fi m ) ■■■ u){x m ,fi m ) 

for every n > 1 and all monotone sequences < x\ < x 2 < • • • < x n , fi- < fi\ < fi 2 < 
■■■ < fi n < fi+ @. 

Suppose that p n ( • ; fix, fi2, ■ ■ ■ , fin) is explicitly known for every distinct fi 
(/i_, fi + ) and that this explicit form can be extended to all fix, t^2, ■ ■ ■ , fJ-n £ K 
words, given an arbitrary n-tuple [fix, fi2, ■ ■ ■ , A*n] G K rf we can produce a p n or, rephrasing 
again, we explicitly know a transformation T : P n — > P n such that 



• • , fin G 
In other 



'M ~W)f = Pn(- -,111,^2, fin)- 



X2) 



It is immaterial whether the definition of u; can be continued to all n 6 
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If, in addition, we know that to is SSC - hence that all zeros of p n reside in (0, oo) - we 
deduce that T maps polynomials with all their zeros in (//_, /a+) into polynomials with all 
their zeros in (0, oo) ||. 

To date, the most efficacious mechanism for the generation of an explicit form of p n has 



been presented in ||], leading to a long list of interesting transformations of the form (1.2). 
Thus, having normalized mo(//) = 1, we asume that there exist real sequences {a n }n£Z+i 
{Pn}neZ+, {7n} ne z+ and {S n } n€Z + such that 



a k 8i-p k 8 t ±Q, k = 0,1,..., £, £ G Z+, (1.3) 



and 



In other words, 



Vr ai + p(3e , 
iHn M = II ; — t-, n G Z + . 

m n+ i(/i) = + f"^ m w (/a), n G Z + . (1.4) 

We call weight functions that possess property (|1.4j ) Mobius quotient functions (MQF). 

It has been proved in []|] that, as long as a; is a Mobius quotient function, the trans- 
formation fll.2j) assumes the form 

i k— 1 n— 1 

x{ e /* n %w n hj(x) \ = e / fc x fc , (i.5) 

fc=0 7=0 i=fc fc=0 



where 



g k (p) = a k + pf3 k , h k (n) = j k + fj,S k , k G 



and where fo, fi, ■ ■ ■ , f n are arbitrary. 

Several MQFs that share strict sign consistency have been presented in || , for example 
uj\{x,p) = x M , ui2(x,[j,) = p x and U2(x,fj,) = x 10 ^^, g G (0,1). Moreover, the theme of 
H was the identification of all possible weight functions a(x) such that u>j(x, p)a(x) is 
MQF and SSC for some j G {1,2,3}. The present paper takes yet another step toward 
our ultimate goal of identifying all MQFs. 

We commence our analysis in Section 2 with the assumption that the coefficients a n , 
Pn, 7 n and 5 n are linear in n. This yields a linear ordinary differential equation with 
polynomial coefficients that must be satisfied by uj and that can be solved explicitly. We 
demonstrate that, up to a transformation of variables, nothing can be gained beyond a 
known transformation from ||. 

In Section 3 we expand our framework to polynomial coefficients a n , (3 n , 7„, 5 n . The 
weight function uj must again obey a linear ordinary differential equation, which we solve 
explicitly in terms of a generalized hyper geometric function. 

Finally, in Section 4 we return to the representation ( |1.5| ), in order to derive an alter- 
native form of the transformation and of the underlying biorthogonal polynomial. Let 



Tlt=i{x-p k ) 
UkZlhk + xSk) 
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We prove that 



\x k , 



p n (x;m,H2, ■ ■ ■ ,/%) = ^F k [qQ,qx,...,q k ]z 

k=0 

where q\ = q(—ctk/Pk), k = 0,1, . . . ,n, and Fq, F\, . . . are a generalization of the familiar 
divided differences. 

2 Linear Mobius quotient functions 

Let us suppose that the moment sequence of the weight function uj obeys (|1.4Q with 

a n = a -na 1 , (3 n = b -nb 1 , 7 n = c -nci, 5 n = d -ndi, n £ Z + . (2.1) 

We further assume that € C 1 [0, oo) and that o;(0, fj,) = — the last assumption 

does not lead to severe loss of generality, since it can be always forced by shifting x i— > x + e 
for some e > 0, smoothly continuing uj for x G [0, e) so that it vanishes at the origin and 
ultimately driving s to zero. We say that a; is a linear MQF.Q 



We rewrite (1.4) in the form 



o 



(in + PnV) x n+ u(x,n)dx = (a n + (3nii) x n uj (x , fl) dx , fl- < /J, < fJL + , (2.2) 
Jo Jo 

multiply by t n and sum for n = 0, 1, . . .. The outcome is the identity 

{(ao — cqx) + n(bo — dox)}oo(x, n)e~ tx dx (2-3) 

= —t {[a\ — c\x) + fi{b\ — d\x)}xuj(x, /j,)e~ tx dx. 
Jo 

Let 

ojt{x,ix) = x l {{ai- C£x) + n(bt- dix)}u(x,fj), £ = 0,1, fi e (fi-,fj, + ). (2.4) 
We can thus rewrite (|2.2| ) as 

C{u (-,n)}(t) = -t£{ui(',n)}(t), 
where £ is the Laplace transform. However, for every / S C , 

£{f}(t) = t£{f}(t)-f(0) 
and we use the smoothness of u> and toi(0, n) = to argue that 

£(i-u 1 (-,A(t) = £{u; (-,fi)}(t). 



\ dx 

Taking inverse transforms, we thus deduce that 



8 

—ui(x,ii) = -uJo(x,u), x > 0, ne{n + ,n-). (2.5) 



The phrase 'linear' refers, of course, to the coefficients a n , /3 n etc., not to oj itself! 
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We finally substitute the explicit values of ujq and u\ from ( |2.3| ), and this yields a linear 
ordinary differential equation for u, 

d 

x{(ai + fih) - (ci + pdi)x} — lo(x, fi) (2.6) 
+ {[(a + pb ) + (ai + /u6i)] - [(c + /ud ) + 2(ci + pd 1 )]x}uj{x,p). 

This can be solved and, in tandem with the initial condition uj(0,p) = 0, results in a 
unique value of u. 

For example, letting 

a n = n, f3 n = l, 7 n = n + 1, (5^ = 1, n£ Z + , 

we obtain 

ao = 0, ai = —1, 6o = 1, b\ = 0, Co = 1, c\ = —1, do = 1, di = 
and ( |2.6P becomes 

-^-v(x,{J,) = - — -u(x,p), 
ox X 

with the solution (up to a multiplicative constant) lj(x,p) = x^ 1 (where x > 0, p > — 1) 
- the Jacobi transformation (with a = 0) from 0. 

With greater generality, we observe that, up to a multiplicative constant, the nonzero 
solution of 

x(a - o x x)y' + (p - p x x)y = 0, y(0) = 0, 



is 



y( x ) = x ~ po/a °(a - a 1 xy o/a °- pi/ai 



- subject, of course, to the above expression being well defined, i.e. to po/&o < and 
either 

Po . Pi . 

> , (To > (J\X 

(To (Tl 

or 

Po , Pi 

< , (To < (J\X 

(TO <Tl 

In our case 

Po = 1 Qq + A^o 
cr ai + pb\ ' 

hence po/ao < is equivalent to 

a + pb 



ai + p6i 

Moreover, 



+ 1< 0. 



Po_ _ Pi _ _ 1 Qq + A^o _ cq + /id 
fj ai + /i^i ci + pd\ ' 
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Particular substitutions readily yield an explicit form of lj. For example, we can use 



" 1 




" 




' 1 




1 ' 


, b = 


, c = 


, d = 







-1 









-1 



and this results in 

Uj(x, fj,) = X f , [I > 1. 

More general is a choice of a, b, c and d so that uj{x,[i) = cx u ^ for some c > and 
a function ^ const. Thus, we need Pq(J\ = pioo and this, after some manipulation, 
results in 

oic = (a -ai)ci, 
Mo = (&o-&i)di, 
(ao&i — ai6o)(aidj — b±ci) = 0. 



Suppose first that ao&i = ai&o- Then 



1 + 



ao 
ai : 



hence f = const, a contradiction. We thus deduce that ai<ii = b\c\. By letting c\ = na\ 
we deduce the explicit expression 



di = biK, c = (ao - aiK d = (b - 
where n ^ is arbitrary - as are ao, a\, bo and &i. We deduce from (2.7) that 



(2.7) 



oi n = 0,0 — ain, 

Pn = b - hn, 

7„ = K(a - ax) - Ka\n, 

S n = K(b - b\) - Khn, 



n £ 



and 



m„(/i) 



„-i 



(ao - airi) + (6 - b\n)p, 
((a - ai) - ain) + ((6 - h) - &m)/x' 

ao + M&o 



ai + /xfei 



We deduce that 

5n(^) 



-(ai + 6i/u) — hn , 

I «i + h/j, 

-refai + &ia) ; h n + 1 

\ ai + bin 



n € 



(2. 



Moreover, aib ^ ao»i readily implies (|L~3l). 
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The function v(n) being strictly monotone, we deduce strict sign consistency of to from 
strict sign consistency of the function x M [12, [HJ. Moreover, 



fc-1 n-l 



( a +b x 
\ ai+bix 



n 9i(x) n hj(x) = (-ir +i K n - k ( ai + b ix r v - a : + z + i , k=o,i,..., n , 

j=0 j=k a\+b\x 

where 

(a) n = a(a+ l)(o + 2) ••• (a + n- 1), a G C, n G Z+, 
is the familiar Pocchammer symbol. Substitution in ( |1.5[ ) thus yields the transformation 

{/ n \ n — k i \ n 

(n«„ + - («o + *») g o (oi+M ;4 +M) } - g 

This, however, can be converted, by a change of variable, to the jacobi transformation 
with a = ||, which can be rendered in the form 

l k =o K + x ) k=0 

The case poai ^ /?i<7o can be resolved explicitly by similar techniques. 



3 Polynomial Mobius quotient functions 

The approach of the last section can be generalized to cater for a more extensive family 
of MQFs. Thus, suppose that the coefficients can be expanded in the form 

oo 

e=o 

oo 

0n = ^h{-n)l, 

1=0 
oo 

In = J2c e {-n) e , 

1=0 

oo 

S n = J2d e (-n) e , 

1=0 

and that the functions 

oo oo oo oo 

a 0) = E = E ^ c ( z ) = E c ^> d ( z ) = E d ^ e 

1=0 1=0 1=0 1=0 

are analytic about the point z = 0. We further assume that uj{ ■ , /i) G C°°[0, oo) and that 



—-^(0,^ = 0, £ = 0,1,.... (3.1) 
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The last condition, in tandem with infinite differentiability, is restrictive. This can be 
alleviated along the lines of the last section, by letting 

J u(x -e,fi) :e<x, 
u e {x, ix) - | 8 ee -(t+i)*- 2 : o < x < e, 

where 9q,Q\,. .. are chosen so as to ensure infinite differentiability at x = e,i and letting 
e — > by the end of the procedure. The situation is simplified if - as we do further in this 
section - it is assumed that a, b, c and d are polynomials, since there it is enough to insist 
that u>( ■ , n) E C r [0, oo) and that (3.1) holds for I = 0, 1, . . . , s — 1, where 

s = maxjdeg a, deg 6, deg c, deg cfj-. 



We again multiply Q2.2| ) by i n and sum up for all integer nonnegative n. This results 
in the identity 

x{c(tx) + [id(tx)}uj(x, n)e~ tx dx = {a(tx) + nb(tx)u(x, fJ,)e~ tx dx. (3.2) 

Jo 

Expanding 

x{c(te) + fid(tx)} and {a(tx) + 
in powers of t, we deduce from (|3.2|) that 



°° />oo 

yV / ^(x,/x)e- te dx = 0, (3.3) 



where u;^ has been given by (|2.4| ), Since (2.4) and (|3.l| ) combine to produce 

^j^(0,/x) = 0, j = 0,1,...,^-!, 



we conclude that 



-tx 



/ / wt(x,n)e tx dx = C^^ut(x,n)j (t), £G 



Taking inverse transforms, (|3.3| ) becomes 

P^(/i) = at + QH/") = Q + ^ 



00 

Y,7^M x ^) = °- ( 3 - 4 ) 



£=0 

Let us set 



therefore 

oj e (x, (Jt) = (Pi{n) - Qe(n)x)x e uj(x, (i), i^lr 



3 The existence of such coefficients follows easily from the Muntz theorem 
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Substitution into (3.4) and differentiation with the Leibnitz rule result in 
Note that both 



1=0 1=0 £=0 



and 



OO f \ OO OO T 

«=0 *=0 ^=0 

are linear in /x. 

We henceforth assume that o, b,c,d S Ps[/i], designating the underlying w as polynomial 
MQF. Therefore ( |3.5| ) is an ordinary differential equation of degree s, equipped with zero 
initial conditions at the origin. 

The equation (3.5) is a special case of the s-degree ODE 

J2(pe-qix)x e yW =0. (3.6) 
1=0 

We express the solution of fl3.6|) in the form 

OO 

y(x) =x 9 ^2 VnX n - 

n=0 

Substitution into ( |3.6| ) and elementary manipulation affirm that, for arbitrary y$, 9 must 
be a zero of 

= (3.7) 

(thus, provided that p s ^ 0, there are s solutions in C, the right number!), whereas the 
coefficients y n are obtained from the recurrence 

., _ SLo(-l)V(-n + l-^ n _ 12 



P(s) := £(-l)^(-z - B) t , Q(x) := ]T(-l)V(-# - x) e . 

i=o t=o 



Let 



Therefore, ( |3.7D is equivalent to -P(O) = 0, whereas ([3^) converts to 

Q(n-1) AQM 1 9 

^^i^ni^. n=l,2,.... 
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Suppose that 



P(x) =p*J[(x -rjj), Q(x) = q*Y[{x - Cj), 

3=1 3=1 

where, without loss of generality, n s = (recall that P(0) = 0). Then 

n s n—1 s 

n p ( £ ) = p* IK 1 - n = ^ n(-o). 

1=1 j=l 1=0 j=l 



therefore 



Vn = — 



n 



j=i\ 



-0) 



ii;; ,n 'io 



-yo, 



n G 



Letting, without loss of generality, yo = 1> we deduce that 



y{x) = x 9 s F. 



s-l 



— (l, — C2> ■ ■ ■ , — Cs! 9* 
1 — 771, 1 — ?72, - - - ,1 -T/a-ll 



a generalized hypergeometric function [14]. 

Of course, we need 9 > s, otherwise the requisite number of derivatives at the origin 



fails to vanish - note that we need only £ = 0, 1, . . . , s — 1 in (p7L 



The last result can be somewhat generalized by observing that the degree of either 
P or Q may be strictly less than s and repeating our analysis without any significant 
changes. 

Theorem 1 Let ct n , (3 n , j n , 5 n G P s [n], n G Z + , and suppose further that cj(-,fi) G 
C s [0,oo) and that d^u(0,fj,)/dx = 0, £ = 0, 1, . . . , s — 1. Then there exist nonnegative 
integers s±, s 2 < s and real numbers 9, Ci) C2, • • • , Csu ^1)^2, • • • , f]s 2 -i^ ^> all functions of 
/i, such that 



u(x,n) = x e Sl F S2 _i 



~(ii —(2, ■ • • ; — Csi! 

1 -7/i,l -772,..., 1 - 77 S2 _i; 



(3.9) 



□ 



We remark that, needless to say, not every Sl F S2 ^i function will do in ( p.9|) , since we 
require io(x,fi) > for all x G [0, 00), /i_ < \i < Moreover, the more interesting case 
is when, in addition, uj is SSC. 

As an illustration of our technique, we consider in detail the case when the quotient 
( |1.4| ) is linear in fj,, i.e. 7 n = 1, 8 n = 0, and assume in addition that s = 2. Therefore 



n G 



«n = oo - 01 7i + a 2 n(n - 1), [3 n = b - bin + b 2 n(n - 1), 
and we have 



P Q (n) = ao + bo/i, 
Pi(/i) = ai + hfi, 
P 2 (fJ.) = a 2 + b 2 /i, 



Qdp) = 1, 

Qi(fi) = 0, 
Q 2 (m) = 0. 
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In other words, 



PoO) 
Pi(m) 

p 2 (m) 



a + ai + ^a 2 ) + (b + 61 + ±6 2 ) /x, 
(ai + a 2 ) + (61 + 6 2 )/U, 
a2 + ^ 



and the differential equation becomes 



i(a 2 + &2/w)s 



+ ((01 + a 2 ) + (61 + 6 2 )/u)x 



1o(m) 

ii(^) 
%(m) 



du>(x, /i) 



1, 
0, 




0, 



().r- 1 ' ' " ' ' /Ac 

+ ((a + «i + ^a 2 ) + (b + 61 + \b<^j fi - x) w(x, /x) 

with zero initial conditions at x = 0. Since 

P(x) = p (/i) + + 0) + ±p 2 (^)(x + 9){x + 6-1) 

and we require -P(O) = 0, we choose 9 = 8(fi) as a root of the quadratic 

|(a 2 + & 2 /i)6> 2 + ((ai + |a 2 ) + (&i + |fe 2 ) + (a + ai + ±a 2 ) + (b + h + |fe 2 ) A* = 0. 

(3.10) 



Therefore 



771 = 1 - 29 - 4 



Pi 

P 2 



and s 2 = 2. On the other hand, Q(x) = 1, therefore s± = 0. Moreover, v = 4/p 2 and we 
deduce from Theorem 1 that 



u(x, n) = x qFi 



— ; Ax 
2# + 4 Pl /p 2 ; ^ 



(3.11) 



Recall that a modified Bessel function can be represented in the form [|T4|] 

1 



I T (z) 



r(r + i) 



r + 1; 4 



Let us assume that /Lt_, a 2 , 6 2 > 0. Then, dispensing with a multiplicative factor, we have 

^,/i) =X^- 2 Pl/P2/ 20+4pi/ p 2 _ 1 



P2 



As long as 9 is real and > 1 we are assured that the initial conditions are satisfied and 
that u>(x, fi) > for all x > 0. 

Let us consider the special choice 

ao = a\ = bo = 6 2 = 0, a 2 > 0, b\ ^ 0. 

In other words, a n = a 2 n(n — 1) and n = —b\n for n S Z + . Letting A = — &i^t/a 2 , ( |3.10| ) 
becomes 

fl 2 + (3 - 4A)# + (2 - 4A) = 0, 
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with the solutions 9 = — 1 and 9 = 4 A — 2. Since we require 9 > 1, we choose the second 
solution and impose A > |. The outcome is 

where /i need be restricted consistently with 61// + §02 < 0. This can be simplified, e.g. 
by considering the variables 




instead of x and u respectively. This results in 

u)(x,p) = x > 0, /l > 0. 



4 An inversion formula for a biorthogonal representation 

The main result underpinning the transformation (|1.5| ) is Lemma 5 of ||, namely that, as 
long as 

n n k—1 n—1 

n (x - ^ = j2 fk n + n ^ + ^,), ^ 

k=l fc=0 j=0 j=fc 

and the moments of u {nx s (/i)} sg ^+ are consistent with ( |l.4| ), the nth biorthogonal poly- 
nomial is 

n 

p n (x;/Ui,U2,...,u n ) = ^2f k x k . (4.2) 

fc=o 

Our purpose in the present section is to express the coefficients fo, fi, ■ ■ ■ , f n explicitly 
in terms of the parameters ui, U2, . . . , fx n . Let 

lJ nU(7 fc +^)' 

Then flO]) yields 

n fe— 1 

$(*) = £/* II "vto- ( 4 - 3 ) 

fc=o j=o 

We let Xj = —a.j/(3j, j = 0,1,..., n and observe that, by (|L~4|), rrij(A^) = for £ = 
0, 1, . . . ,j — 1. Therefore, letting x = Xe in ( |4.3D , we obtain 

t k-i 

^ = E/fcIT m i(^)' * = 0,l,...,n, (4.4) 
k=o j=0 

where q~i = q(Xg), £ = 0, 1, . . . , n. 
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We assume that rrij(A^) ^ for £ > j. Letting I = 0, 1, . . . in (|Q|), we obtain 



fo 
fi 

h 



go, 

gi - go 

m (Ai)' 

g~2 - go _ gi - go 
m (A 2 ) mo(Ai) 



mi(A 2 ) 

etc. To obtain the general form of f k , we need to generalize the concept of divided 
differences. Given a sequence {re}t£Z+, we define 

^oM = ri, 

Fk-i[rt, ■ ■ ■ ,r e+k ^ 2 ,re+k - F k _i[r e , . . . ,r i+k -i] 



F k [r e ,r e+1 , . . .,r e+k ] 



Lemma 2 f k = F k [q Q , q\, . . . , q k ], k = 0, 1, 



Wfc-l(A^+fe) 



fc = l,2. 



,n. 



Proof. By induction on The statement is certainly true for £ = and we assume it 
for £ — 1. Therefore, it follows from (4.4) that 

i-\ k-1 t-l 

Qe = go + E F k[qo, • • • , gfc] II m ( A i) + d t II m i( A ^)- 

k=l j=0 j=0 



Therefore, 



E 



, F fc [% . . . , q k ] 



Fi[qo,< 



E 



* -Ffc [go • • • , gfc] 



We will prove by induction on s = 1,2,.. 

^s[g0i • • • 3 Qs-i, 



2 that 



ft 



nJ-=>i(A<) 



E 



-v -Ffc[g0,gl, •••,<&] 



I&WA«) 



(4.5) 



This is certainly true for s = 1 and, assuming that (4.5) is valid for s > 1, we have 



F s [(fc, • • • , gs-l,%] - F s [q , ... , g s 



nt>i(A<) 



i-l 

E 



^fc[g0,gl, ...,%] 



k=s+l Hj 



^+l[g0, • • • ,gs,gs+l] _ ^4 /"/,- [(/(). gi • • • • 



n 



This proves that (4J>) is true. We now let s = £ — 2, therefore 

Fi-i[qo, • • • , qe-2, Qe] — -F)-i[go, ■ • • , qe-2, Qe-i] 



ft 



m^_i(A^) 



Fe[q ,Qi, 
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This completes the proof. 



□ 



Theorem 3 Let the moments of u satisfy ([O]) an d assume that (3% ^ 0, £ £ Z + , and 

Oijfti aiPj for all £,j € Z + , I ^ j. Set 





Then 



p n (x; fix, (i 2 , ■ ■ ■ ,Mn) = ^-Fjfc[go,gi,---,9fc]a; fc - 

fc=0 



(4.6) 



Proof. Follows at once from the lemma, since the above definition of {qe}£ =0 is consis- 
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